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We calculate, at the one-loop level, penguin diagrams for improved staggered fermion operators 
constructed using various fat links. The main result is that diagonal mixing coefficients with penguin 
operators are identical between the unimproved operators and the improved operators using such 
fat links as Fat7, Fat7+Lepage, Fat 7, HYP (I) and HYP (11). In addition, it turns out that the 
off-diagonal mixing vanishes for those constructed using fat links of Fat7, Fat7 and HYP (II). This 
is a consequence of the fact that the improvement by various fat links changes only the mixing 
with higher dimension operators and off-diagonal operators. The results of this paper, combined 
with those for current-current diagrams, provide complete matching at the one-loop level with all 
corrections of 0(g 2 ) included. 

PACS numbers: 11.15.Ha, 12.38.Gc, 12.38.Aw 



I. INTRODUCTION 

Decoupling of heavy particles in the standard model 
leads to the low-energy effective Hamiltonian through 
OPE (operator product expansion) which achieves a fac- 
torization of short- and long-distance contributions. The 
effective Hamiltonian of our interest consists of AS = 1 
four-fermion operators which contain information at low 
energy and the corresponding Wilson coefficients includ- 
ing the short-distance physics at high energy. The low en- 
ergy effects of the electroweak and strong interactions can 
be expressed as matrix elements of the four fermion oper- 
ators between hadron states. Since the matrix elements 
involve non-perturbative QCD physics at low energy, lat- 
tice QCD is well suited for their calculation. The main 
advantage of using lattice QCD is that it provide a first 
principle, non-perturbative estimate. Different fermion 
discretizations such as Wilson, staggered, domain-wall 
(DW) and overlap have been used to study the matrix 
elements and to control their statistical and systematic 
errors. 

In this work, we adopt staggered fermions to explore 
the renormalization property of the operators relevant 
to calculating the weak matrix elements. Staggered 
fermions are an attractive choice for the calculation of 
weak matrix elements. They are computationally effi- 
cient and so simulations with three dynamical flavors are 
already possible with relatively light quark masses using 
the AsqTad staggered fermions 0] . They preserve enough 
chiral symmetry to protect operators of physical interest 
from mixing with others of wrong chirality. By construc- 
tion they retain four tastes of doublers for each lattice 
field, which in itself is not a problem. Their main draw- 
back is that the tastes symmetry is broken at non-zero 
lattice spacing and is restored in the continuum limit. 
At non-zero lattice spacing, quark-gluon interactions vi- 
olate the taste symmetry. This has two consequences: 
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first, there are large 0(a 2 ) discretization errors in hadron 
spectrum and weak matrix elements and second, some 
one-loop corrections are so large that matching factors 
differ significantly from their tree-level value of unity. 

Both problems are manifestly alleviated by improving 
staggered fermions using "fat" links 0, • Taste sym- 
metry violations in the pion multiplet are substantially 
reduced 0, 0| and one-loop corrections to the match- 
ing factors for four-fermion operators are reduced from 
as large as 100% down to ~ 10% Based upon the 
analysis of Ref . 0, , we figured out that the perturba- 
tive corrections is smallest for a particular type of fatten- 
ing, "HYP" and "Fat 7" smeared links. Hence, we choose 
the "HYP/Fat7" improvement scheme for our numerical 
study on the weak matrix elements @- 

Recently, there have been elaborate efforts to do the 
higher loop calculation for highly improved staggered 
fermions using the Luscher and Weisz method [13. They 
automated the generation of the Feynman rules for essen- 
tially arbitrary complicated lattice actions such as the 
AsqTad type [TJ and some even developed an algorithm 
which generates Feynman diagrams automatically 12]. 

An essential step for using lattice QCD is to obtain 
the relationship between the continuum and lattice four- 
fermion operators. There are two classes of Feynman 
diagrams at the one-loop level: (1) current-current dia- 
grams and (2) penguin diagrams [ljjllj]- At the one-loop 
level, the contribution from the current-current diagrams 
and the penguin diagrams can be treated separately. In 
the case of the current-current diagrams, the matching 
coefficients for the operators constructed using improved 
staggered fermions (HYP/Fat7) were presented in 

Here, we focus on penguin diagrams in which one of the 
quarks in the four-fermion operator is contracted with 
one of the anti-quarks to form a closed loop. The main 
goal of this paper is to calculate the penguin diagrams 
for improved staggered fermion operators constructed us- 
ing various fat links and to present the matching formula 
between the continuum and lattice operators at the one- 
loop level. The results are compared with those for the 
unimproved staggered fermions given in |14| . The role of 
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improvement is reviewed from the standpoint of operator 
mixing. The main results of this paper, combined with 
those of current-current diagrams given in Ref. 0, pro- 
vide a complete set of matching formula for e'/e at the 
one-loop level with all the g 2 corrections included. 

This paper is organized as follows. In Sec. [HJ we de- 
scribe our notation for the improved staggered fermion 
operators and various fat links. We also explain the Feyn- 
man rules. In Sec. II I II we explicitly calculate penguin 
diagrams step by step and the main results are summa- 
rized in a theorem. We close with some conclusions in 
Sec. El 

A preliminary result of this paper appeared in |l5j . 



II. NOTATION AND FEYNMAN RULES 

Basically, the improved staggered fermion action has 
the same form as for the unimproved staggered fermions 
[Tfil | , except in that the original thin links £/ M are replaced 
by fat links V^. 

S = a 4 J! ^~^2 r l^(ri)(x(n)V^(n)x{n + fl) - 



X(n + p)Vhn)x{n) \ + mx(n)x(n) 



(1) 



where n = (m, 112, 713, 114) is the lattice coordinate and 
Vn( n ) = ( — ^"f- 1 . Here, the fat link repre- 
sents collectively Fat7 [I]}, Fat 7+Lep age [11. HYP links 
@ and SU(3) projected Fat7 (Fat7) links[|. Th e de- 
tailed definitions of various fat links are given in the 
original references and so we do not repeat them here. 
In addition, MILC and HPQCD collaboration developed 
and have used the AsqTad staggered action 0, 0, fill > 
whose gauge action is a one-loop Symanzik improved ac- 
tion and whose fermion action contains a Naik term in 
order to remove the 0(a 2 ) discretization errors. Out of 
various choices, HYP and Fat7 suppress, in particular, 
the taste changing gluon interactions efficiently and re- 
duce the taste symmetry breaking in the pion spectrum 
significantly 0, compared with others. In addition, it 
turns out that the one-loop corrections are smallest for 
HYP and Fat7 compared with others and that they 
possess several nice properties in renormalization which 
are explained in . If the goal of the improvement were 
to minimize the 0(a 2 ) discretization error and to achieve 
a better scaling behavior through the Symanzik improve- 
ment program, it would be natural to choose the action of 
the AsqTad type and improve the operators correspond- 
ingly. However, our improvement goal is to minimize the 
perturbative correction as much as possible (if possible, 
down to less than 10% at the one loop level), which, in 
fact, turned out to be the same as to minimize the taste 
symmetry breaking effect |l9j . Hence, we have chosen 
HYP/Fat7 improvement scheme for our numerical study 
on e'/e mainly because it serves better to our purpose of 



improvement. In this paper, we adopt the same notation 
for fat links as in 00, 

In order to construct the four spin component Dirac 
field, we adopt the coordinate space method suggested in 
|20j. In this method, we interpret 16 staggered fermion 
fields (x) of each hypercube as 4 Dirac spin and 4 fla- 
vor (=taste) components. The continuum limit of the 
staggered fermion action on the lattice corresponds to 
QCD with four degenerate flavors [Nf = 4) [|(j. There 
are numerous choices to transcribe the lattice operator 
for a given continuum operator pH I22I I23L We adopt 
the same convention and notation as in |24| except for 
the gauge links. We denote the gauge invariant bilinear 
operators as 

[SxF] = -r^X[xa(jM) (7s®£f)ab XbiVBt)] 



A,B 



V ab (y A ,y B ) 



(2) 



where S and F represent spin and flavor (=taste) respec- 
tively and correspond to one of the following; S (scalar), 
(vector), T^ v (tensor), (axial) and P (pseudo- 
scalar). Here, V(yA,ys) is a product of gauge links that 
makes the bilinear operator gauge invariant. The link 
matrices V(yA,yB) are constructed by averaging over all 
of the shortest paths between tja and ys , such that the 
operator [75 x £p] is as symmetric as possible: 



V{y A ,yB) = ^^2v{y A ,y A + A P1 ) - ■ ■ 
p 

V(y A + A P1 + A P2 + A P3 , y B ) , (3) 

where P is an element of the permutation group (1234) 
and 



A,, 



(4) 



For the four-fermion operators, we use the same notation 
as the bilinears but need to distinguish between color one 
trace and color two trace operators. 

1 



[S x F][S' x F']j ee -j 



N 2 ^ 

1 A,B,C,D 



[S x F][S' x F'\n 



\x{Va) {is®£f)ab x{Vb)] 
■[x{Vc) hs' ®Cf')cd x(vd)] 
■V(yA,VD) V(y c ,y B ) (5) 

N 2 ^ 

f A,B,C,D 



Mva) (7s ® £f)ab x(yB)} 
■ixivc) (is* ®Cf0cd x(vd)] 
■V{va-,Vb) V{yc,VD) (6) 

Here, note that the sub-indices /, II represent the color 
one trace and color two trace operators respectively. 
There are two completely independent methods to con- 
struct operators on the lattice using a Fierz transforma- 
tion: one spin trace formalism and two spin trace for- 
malism pi} . In this paper, we choose two spin trace 
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formalism to construct the lattice operators and it is also 
adopted for our numerical study on e'/e. 

The Feynman rules for the unimproved staggered 
fermions are standard and presented in (Til l2ll Y2M. 12M 
l26l I27I Esl ] . Here, we use the same notation for Feynman 
rules as in [2l], |24j and we do not repeat them. By intro- 
ducing fat links, we need to change the Feynman rules. 
These changes in the Feynman rules originating from fat 
links are given in [|| Q • Here we adopt the same nota- 
tion and Feynman rules as in 0, Q • We explain only the 
essential ingredients for a one-loop calculation. 

We define the gauge fields of the thin link and fat links 

as 

U^x) = cxp (iaA^x + (7) 

V^(x) = cxp (iaB^x + -ft)) (8) 

Here, note that corresponds to various fat links. We 
call a "fat gauge field" and a "thin gauge field" . 
This fat gauge field can be expressed in terms of thin 
gauge fields as follows: 

oc 
ra=l 

= BW + B( 2 '+0(i 3 ) 

Here, represents a term of order A n . Theorems 1 

and 2 of say that the linear term is invariant under 
SU(3) projection and that since the quadratic term is 
anti-symmetric in thin gauge fields, its contribution van- 
ishes at the one-loop level [33. Hence, at one loop level, 
the renormalization of the gauge invariant staggered com- 
posite operators constructed using fat links of HYP type 
can be done by simply replacing the propagator of A^ by 

that of B^ . This simplicity is extensively used to cal- 
culate the one-loop correction to the improved staggered 
operators 0,0- 

The linear term B^ can be expressed in momentum 
space as 

The details of the blocking transformation for the fat 
links are contained in h fJlV (k), which is given in 0,0 for 
various fat links. 

K„{k) = S^D^k) + (1 - 5^)G MV (k) 
D^k) = l- dl J2t + d2 £ &l 

-d 3 sls 2 p sl -d 4 ^2st 

G^(k) = s^G.Ak) (9) 
G^(fc) = dl -d 2 (s2p+ 2 Sl) +d 3 S -^+d 4 sl (10) 



Here, the coefficients di distinguish the different choices 
of fat links. 

1. Unimproved (naive): 

di = 0, d 2 = 0, d 3 = 0, d 4 = 0. (11) 

2. Fat 7 links: 

di = 1, d 2 = 1, d 3 = 1, di = 0. (12) 

3. HYP (I) links: 

d l = (2/3)01(1 + 02(1 + 03)), 

d 2 = (4/3)oi0 2 (l + 2o 3 ), 

d 3 = 8010203, 

d 4 = 0. (13) 

We consider two choices for the o^. The first was 
determined in |5j using a non-perturbative opti- 
mization procedure: oi = 0.75, a 2 = 0.6 03 = 0.3 
(we call this choice "HYP (I)"). This gives 

rfi=0.89, d 2 =0.96, d 3 = 1.08, d 4 = 0. (14) 

4. HYP (II) links: The second is chosen so as to re- 
move 0(a 2 ) flavor-symmetry breaking couplings at 
tree level. This choice, oi = 7/8, a 2 = 4/7 and 
03 = 1/4 (we call this choice "HYP (II)"), gives 

dl = 1, d 2 = 1, d 3 = 1, di = 0, (15) 

i.e. the same as for Fat7 links. 

5. Fat7+Lepage (0(a 2 ) improved links): 

di =0, d 2 = 1, d 3 = 1, di = 1, (16) 

For later convenience, we name the SU(3) projected Fat7 
scheme "Fat 7" . 

In |§] , we studied various fat links from the standpoint 
of renormalization to improve staggered fermions and it 
turns out that at the one-loop level, the renormalization 
of staggered fermion operators is identical between Fat 7 
and HYP (II). In addition, it was explained that SU(3) 
projection plays a role in tadpole improvement for the 
staggered fermion doublers. 

III. PENGUIN DIAGRAMS 

Here, we study penguin diagrams in which one of the 
quarks in the four-fermion operator is contracted with 
one of the anti-quarks to form a closed loop. The main 
goal is to calculate the penguin diagrams for improved 
staggered fermion operators and to provide a matching 
formula between the continuum and lattice operators at 
the one-loop level. 
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FIG. 1: Penguin diagrams for the staggered fermions 



On the lattice, the gauge non-invariant four-fermion 
operators such as Landau gauge operators mix with 
lower dimension operators which are gauge non-invariant 
[Tij . It is required to subtract these contributions non- 
perturbatively. However, it is significantly harder to ex- 
tract the divergent mixing coefficients in a completely 
non-perturbative way. This make it impractical to use 
gauge non-invariant operators for the numerical study of 
the CP violation. Hence, it is necessary to use gauge- 
invariant operators in order to avoid unwanted mixing 
with lower dimension operators. For this reason, we 
choose gauge-invariant operators in our numerical study. 

In the staggered fermion formalism there are four pen- 
guin diagrams at the one-loop level as shown in Fig. ^ 
These diagrams of penguin type can mix with lower di- 
mension operators in addition to four-fermion operators 
of the same dimension or higher. The mixing coeffi- 
cients with lower dimension operators are proportional 
to inverse powers of the lattice spacing. The perturba- 
tion, however, is not reliable with divergent coefficients. 
Hence, we must use a non-perturbative method to de- 
termine them and subtract away the lower dimension 
operators. In the case of mixing with operators of the 
same dimension, the perturbativc calculation is expected 
to be reliable as long as the size of the one-loop correc- 
tion is small enough, which can be achieved naturally by 
improving the staggered operators using fat links. 

In Fig. n diagrams (a) and (b) have their correspon- 
dence in the continuum and diagrams (c) and (d) do not 






(e) (f) 

FIG. 2: Diagram identity: (c) + (d) = (e)+(f). 



have any continuum correspondence. However, diagrams 
(c) and (d) play an essential role in keeping the gauge 
invariance in the final sum. In other words, the gauge 
invariance is broken without them. 

First, we overview the role of each diagram in the 
gauge invariance and present the details later. Basically, 
the contribution from diagrams (c) and (d) can be re- 
expressed as a sum of two separate terms: diagram (e) 
and (f) in Fig. |2 We observe that the sum of diagram 
(a) of Fig. ^and diagram (e) of Fig. [21 produces bilinear 
operators in a gauge invariant form as shown in Fig. [3J It 
turns out that the contribution from diagram (b) in Fig.^ 
and diagram (f) in Fig. [2] leads to four-fermion opera- 
tors of our interest in a gauge invariant form, which are 
typically called "Penguin diagrams" in the literature as 
shown in Fig.^ Once more, we emphasize that diagrams 
(c) and (d) (or equivalently (e) and (f)) are essential to 
keep the gauge invariance. 



A. Penguin operators on the lattice 

In order to construct a lattice version of continuum 
penguin operators, we need some guide lines, because the 
staggered fermions carry four degenerate tastes by con- 
struction, unlike the continuum fermions. Hence, a closed 
loop of staggered quarks contains four degenerate tastes 
running around it, rather than a single quark, which 
needs to be normalized properly by 1/Nf — 1/4. Pen- 
guin diagrams occur for operators which belongs to octet 
irrep of the continuum flavor SU(3) symmetry. Apart 
from the overall factor, the calculation is identical for all 
the penguin operators of our interest 0,H3- Therefore, 
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(b) + (f) 
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FIG. 3: Bilinear mixing. 



we may choose the following operator as a representative 
without loss of generality: 



U S',S 



-jLatt 



where the operators are defined as 



(17) 



O^ 



1 



s?f>,sf = ir^iXshs* 



1 €F')Xd][Xu("/S ® £f)Xu] 



Here, we adopt the two spin trace formalism [2l(, color 
contractions with gauge links are dropped for brevity, 
and the subscripts s,d,u represent the continuum quark 
flavors. Here we select the continuum flavors so that 
there is only one possibility of up quark contraction to 
form a closed loop, leading to a penguin diagram. The 
bilinear with strange and down quarks behaves as a spec- 
tator in the calculation. This choice of continuum flavor 
assignment was suggested originally in [l4|. 

Penguin diagrams of the above operators lead to mix- 
ing with the same class of SU(3) octet operators (penguin 
operators): 



(18) 



q 

^2 iXshs' ® ^F')Xd]J2iXq(ls ® Zf)x 9 ] (19) 



U S'F'SF 



N f 

J q 

where the sums run over the active light flavors such as 
u, d, s. 



(i) <M 

FIG. 4: Mixing with penguin operators 



We adopt the same notation as in 0, [24| to incorpo- 
rate two different color contractions: 



O = 



On 



(20) 



The contribution from the color two trace operator, On 
vanishes in the penguin diagram because it is propor- 
tional to Tr(T a ) = (T a is the SU(3) group generator 
which is traceless). Hence, it is sufficient to work only 
on the color one trace operator, Oj. This simplifies the 
matching formula. 



(0 Con t)/ = (0 Lat t)/ + £ ^ ^ . 



vLatt \ 



(21) 



where i, j includes both the spin and taste indices, and 
the projection vector is 



1 2' 



(22) 



The Zij have two separate contributions: one from the 
continuum operators and the other from the lattice op- 
erators. 



iy r^Cont r^Latt 



(23) 



The main goal of this paper is to calculate , since 
Zf° ut is well known jajj. 
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B. Feynman diagrams 

Using the Feynman rules described in Sec. [HI we ex- 
press the analytical form of the Feynman diagrams pre- 



sented in Fig. ^ Diagram (a) of Fig. ^ contains only 
mixing with bilinears and combined with diagram (e) of 
Fig. [21 makes the bilinears gauge invariant as shown in 
Fig. El 



G(a) = -^rK (a) ( 24 ) 



N f - " 



Hl {a) = S cd (75' ® Hf>) c , d , ^ J2 (rrs ® i F ) AB ■ lf B (25) 

AB 

I { f B = [ exp(ip-(A-B))lS F (p)] BA (26) 
Jp 



Diagram (b) of Fig. combined with diagram (f), Fig. 01 Using the Feynman rules described in Sec. ITTl we 
leads to mixing with four-fermion operators as shown in can express the diagram (b) as follows: 



G {b) = ~ J H*>V{k)'D%(k)-Vj(p + Trc,-q-*D,-k) (27) 

Hj>W(k) = (-ig) Ti d (7ST'®^wiE^®^ B exp (i| • (A + B)) ■ lf B {k) (28) 

AB Z 

I A b) B (k) = / exp(ip-(A-Bj) {[S P (p-±)]-coB(pjJj^.[S F (p+±)]} BA (29) 



Diagram (c) and (d) of Fig. ^ contains mixing with bi- These two diagrams play an essential role in keeping the 
linears and four-fermion operators as shown in Fig. |21 gauge invariance of the final results. 



G (c) = -j^-H^ik) (30) 



H^(k) = ( +i 0)T/ d i£( 7s ®£F) * Y.^~ l ) C ' A '^®^) A , B ,{-l) D - B \B-A'),f% B {k)lf B (31) 



16 ~* ^ ' A *16 „ 

A,B A'.B 



G 



(d) 



-±-HW{k) (32) 
H^ d) (k) = (+ig) T/ rf l ® &)ab^ E (-!) C ' A '(7S' ® WM^V - ^ (33) 



16^ v '" ~~ '^16 , 



C. Gauge invariance 



As shown in Fig. |3 diagrams (c) and (d) of Fig. 

can be expressed as a sum of diagrams (e) and (f). As a 
consequence of this, diagrams (b) and (f) lead to a gauge- 



invariant form of penguin operator and diagrams (a) and 
(e) lead to a gauge-invariant bilinear operator. 



First, let us derive the relationship shown in Fig. |2] 
We begin with a simple identity for diagram (b): 
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^2sin(|).^)(fc) 



1 



AS 



exp(ifc ■ A) — exp(ifc • B) 



r(«) 
MB 



(34) 



r 



The continuum correspondence of the above Eq. (|34|) van- 
ishes but the RHS of Eq. I|34[l is not zero. Then, what 
is going on with the gauge invariance? It turns out that 



term, which is the focal point of this section. 

The same kind of identities for diagrams (c) and (d) 



there exists an additional term which cancels out this can be expressed in terms of diagrams (e) and (f). 



J2, 2 Sm (^).fe< c >(fc)+ff^ (35) 



Here, the H^ e '(k) and H^'ik) are defined as 



H'/ e) (k) = (+ig)T! d ± Y,(- 1 f'' A '^^^) A , B/ -(-l) D '- B \B'-A')^, BI (k)-^Y,^^MABl 



A'B> 



A.B 



H^f\k) = (+ig) Tl d ( 75 , ® £ F ,) C , D , -^£(7s® Cf) ab {B - A),f AB (k) ■ J« 



A,B 



a) 
AB 



(36) 
(37) 



r 



A derivation of the identity given in Eq. I|35|l is given in 
Appendix El 

Note that Hl' [e \k) can be factorized into two bilin- 
ears: one bilinear has a single gluon emitted and the 
other bilinear forms a closed fermion loop identical to 
that of diagram (a). Unlike the continuum operators, the 
staggered operators are non-local and gauge links must 
be inserted between the quark and anti-quark fields to 
make them gauge invariant. Thus, at O(g), we must 
have H^ e \k) to keep the gauge invariance of the spec- 



tator bilinear. Of course, this is completely a lattice arti- 
fact, which will vanish in the limit of zero lattice spacing, 
a = 0. Let us explain this bilinear mixing in detail in the 
next subsection. 



D. Bilinear mixing 

Here, we combine diagrams (a) and (e) and show that 
they form a gauge-invariant bilinear. 



H i,(a) + H Ue) {k) = -J2 (-1) C '' A '(75' ® Zf') a , bi {-1) D '- B ' ■ \s cd + (iga) T^B 1 - A%f% B ,{k) 



A',B' 



16 



(38) 



A,B 



Here, the term enclosed in square brackets is nothing the gauge coupling g. Note that the closed loop part 
but an expansion of gauge link V{da,Db) in powers of 
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behaves as a constant defined as 



X (a+e) = ^ H (7S ® fr) 



r(«) 
AB 1 AB 



(39) 



A,B 



Hence, the combined result of diagrams (a) and (e) be- 
comes 



X 



(a+e) 



x{va){is> ®ZF>) AB V{y A ,yB)x{yB) (40) 



Naturally, the next question would be what are the 
possible spins and flavors for (7s<8>£f) in the closed loop. 
Since we are talking about the vacuum diagram, it must 
have the vacuum quantum number. Hence, the natural 
choice would be that (75 <g> £p) = (1 <g) 1). 

Another question would be how reliable is the coef- 
ficient X( a + e ). Since we are talking about the diver- 
gent coefficient, hence, the perturbative determination of 
X( a+e ) is highly unreliable because the contribution from 
the truncated terms are also divergent as we approach to 
zero lattice spacing. Hence, this coefficient must be de- 
termined using a non-perturbative method. 



level. First, we address the issue of the gauge invariance. 
Basically, we want to show how the RHS of Eq. I|34|) 

i(t) 

cancels out. From the definition of ■ ' (k) given in 
Eq. i|37|) it is easy to show the following Ward identity: 



^2sm(f).(^)(fc) + J ffW)(fc)) =0 (41) 



This illustrates that, unlike the continuum, we need di- 
agram (f) to keep the gauge invariance of diagram (b), 
which is a pure lattice artifact originating from using 
staggered fermions when constructing operators on the 
lattice. 

Now, we turn to explicit calculation of the mixing with 
penguin operators. First, we define I^L in terms of I A a B 



in Eq. as follows: 



Ac) 

AB 



r(«0 



/ exp(ip ■ (A - B)) 
Jp 

■{[S F (p)]'C08(p^^[S F (p)]\ ,( 42 ) 

l. J BA 



E. Penguin operator mixing 

Here, we want to present the main results of Penguin 
diagrams: mixing with penguin operators at the one- loop 



Using I A %, we can collect diagrams (b) and (f) into such 
a form that all of the nice features necessary for the 
gauge invariance and the infra-red behavior are visible 
at a glance. 



= H.9) TU ( 7 S'®6-W Jq E -hsSWAB exp (i~ ■ (A + B)) ■ (l A % (fc) - ) (43) 



This result is identical to that originally presented in [T^ . 
Regarding the infrared behavior, note that 



l^jAB(k) = I. 



k^O 



AB 



Hence, this confirms that i?^ + (&) vanishes when k = 
0. In addition, the quark mass behaves as an infra-red 
regulator such that the integrals of I AB (k) and I AB are 



infra-red safe. As shown in Eq. 1)41(1. H) 



I,{b+f) 



(fc) also 



satisfies the Ward identity coming from gauge invariance: 



;sin(^).ff^+/>(fc) 



(44) 



for all values of k. 

Now we turn to explicit calculation of Hp[^ b+ ^\k) in 
such a form as we can use for matching to the continuum. 



The mixing contributions comes not only from k ~ but 
also, in principle, from k ~ ixja ■ A for any arbitrary 
hypercubic vector A. Here, we first consider k ~ as in 
the continuum. 

n^ s \k) = (-ig)T^^7W^cD 

~J2^ s ®^ab[P^)]ba (45) 



[P»(k)]BA 



A,B 

i 



im 
'(4^)2 

I 2 

' (4tt) 2 



a 

' y^feafc 2 - k^k a ){^ a ® l) BA • I c 

ex 

(46) 
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where the lattice-regularized finite integrals I a , It and I c 
are given in Appendix [BJ Note that this result for the 
improved staggered fermions is identical to that for the 
unimproved staggered fermions presented in [13]. This 
equivalence will be discussed in detail later when we 
present Theorem 1. 

The first term in Eq. 1|46[) describes mixing of the four 
fermion operator with a bilinear with gluon emission: 



\Xs (7S' ® Xd][Xu {Ovp ® iv$) Xu] 



[Xs (75' ® 6 



F vp Xd] 



(47) 



This corresponds to mixing with a dimension 5 opera- 
tor. From the standpoint of physics, this mixing belongs 
to a class of unphysical operators because none of the 
operators of our interest possesses a flavor structure of 

Similarly, the second term in Eq. l|4t)|) corresponds to 
mixing of the four-fermion operators with a bilinear: 



\Xs {lS> ® Xd][Xu {(J pa ®l)Xn 

— ► m[Xs (is' xd] 



(48) 



This represents mixing with a dimension 6 bilinear which 
has its correspondence in the continuum. However, this 
operator vanishes in the chiral limit and so it corresponds 
to a higher order in the chiral perturbation. In addition, 
the spin structure of the tensor does not appear in the 
original set of operators of our interest and so this mixing 
can occur, if any, at order g . Hence, by convention, this 
term is dropped from the analysis [lj, |23 . 

It is the third term in Eq. (|46() that corresponds to 
mixing with penguin operators. To complete the opera- 
tor construction, the gluon needs to be connected to the 
external fermion line as follows. 



G 



(b+f) 



Vj{p + ir c ,-q-n D ,-k) (49) 



where Vj corresponds to the fermion vertex emitting one 
gluon in |2l|. Here, Dj^l(k) represents the gluon prop- 
agator of the thin or fat links which can be collectively 
written in terms of fc M = 2sin(fc M /2) as 



a,0 



k 2 



(1-A) 



k a kp 

\k 2 ] 2 



(50) 



in a general covariant gauge, where k 2 = kp. It turns 
out that only the diagonal part of hp a (k) contributes 
mainly because the off-diagonal term of h^a^k) is pro- 
portional to kpk a — > kpk a and so the contribution from 
the off-diagonal term of hp a {k) vanishes due to a simple 
identity: (Sp a k 2 — kpk a ) ■ k^ = 0. 

($fii/k k fj,k v )h fj, a {k) = {Spuk k^k^S pahfi^ik) (51) 



In addition, since /i MM (fc) = 1 + 0(a k ), the same iden- 
tity also guarantees that the gauge fixing term propor- 
tional to (1 — A) also vanishes for the leading term in the 
limit of k — > 0, which insures gauge invariance. In sum- 
mary, we can claim that in the low momentum limit of 
k -> 0, 



{Sp a k 2 - kpk a ) ■ D"{k) = 5 u 5^[h^{k)] 2 



1 



'(fifiak kpka) (52) 

Another important ingredient is that the k^k a part of 
{5p a k 2 — kpka) cannot contribute in the on-shell limit. 
Using the equations of motion for the staggered fermions, 
we can prove that regardless of quark mass, 



J2sm(^-)-Vj(p,-q,-k) = 



(53) 



Hence, in the limit of small momentum, the S^k^ka term 
vanishes by the equations of motion. In addition, in the 
low momentum limit, 



exp (i-j ■ (A + B)j = 1 + 0{kc 



(54) 



and we are interested only in the leading term, since the 
contribution from the remaining higher dimension oper- 
ators is supposed to vanish as we approach to the con- 
tinuum (a = 0). 

Now we can simplify f<\ defined in Eq. I|49|l at small 

k, 



CD 



Ss,ft 8 F ,i [hpp,(k)} 



(55) 



where k = q — p is strictly on shell. From the above 
Eq. I)55|). we can derive an interesting theorem: 

Theorem 1 (Equivalence) 

At the one-loop level, the diagonal mixing coefficients 
of penguin diagrams are identical between (a) the unim- 
proved (naive) staggered operators constructed using the 
thin links and (b) the improved staggered operators con- 
structed using the fat links such as HYP (I), HYP (II), 
Fatl, Fatl+Lepagc, and Fat7.pl 

Proof 1.1 In the case of the unimproved staggered oper- 
ators, hpp(k) — 1 by definition. The improvement using 
the fat links such as HYP, Fat 7, Lepage+Fat7, Fat7, in 
general, leads to as defined in Eq. (|10fl . The role of 

the additional terms proportional to d t is to suppress the 
high momentum gluon interactions at the cut-off scale 
(it /a). Hence, by construction, these additional terms 
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cannot change the dispersion in the low momentum re- 
gion. In other words, in the limit of k — > 0, 

= 1 + 0(k 2 a 2 ) (56) 

Here, the 0{k 2 a 2 ) term corresponds to higher dimen- 
sion operators, whose contribution vanishes in the limit of 
a = 0. Hence, this term is irrelevant to the penguin mix- 
ing of our interest. It is only the leading term of Eq. ljST)|) 
that contributes to the mixing with penguin operators. 
The [/i MM (/c)] 2 term is, if any, the only possible source of 
difference introduced by the improvement. However, the 
contribution from [h^^k)} 2 is, by construction, identical 
in the low momentum limit before and after the improve- 
ment. Therefore, this leads us to the conclusion that the 
mixing coefficients with penguin operators must be iden- 
tical for the staggered operators constructed using both 
thin links and fat links such as HYP, Fat7, Lepage+Fat7, 
and Fat7. This completes the proof of the theorem. 

As a consequence of Theorem 1, it is trivial to obtain 
the diagonal mixing coefficients from Eq. I|55|) . 

= ~Ic (57) 

where i represents (7^ <g) 1). This is our final result. 



F. High momentum gluons and off-diagonal mixing 



when G p ^ v = 1 in at least one transverse direction. The 
vertex also vanishes when C v — \. Hence, this off- 
diagonal mixing is absent at the one-loop level. Since 
we adopt either Fat7 or HYP (II) as our improvement 
scheme in our numerical study, we rejoice in this absence 
of unphysical off-diagonal mixing when we analyze the 
data. 

In the case of the improvement using HYP (I) and 
Fat7+Lepage, h^(k + He) does not vanish exactly but 
it is significantly suppressed. Correspondingly, the off- 
diagonal mixing is similarly suppressed. 



G. Tadpole improvement 

In [24| , the procedure of tadpole improvement for the 
staggered four-fermion operators is presented. The tad- 
pole improvement factor is given, basically, in powers of 
uq. In perturbation, the contribution from this is of order 
g 2 so that only the tadpole improvement of the original 
operator at the tree level can contribute at the g 2 or- 
der. Hence, the one-loop result for penguin operators 
are of order g 2 and the tadpole improvement can change 
the result only at the order of g 4 . In other words, the 
tadpole improvement corrections included when we cal- 
culate the current-current diagrams are complete at g 2 
order and so there is no additional correction from the 
tadpole improvement to the penguin diagrams. This ar- 
gument holds valid not only for the diagonal mixing but 
also for the off-diagonal mixing. 



By construction, gluons carrying a momentum close 
to k ~ ix J a are physical in staggered fcrmions and lead 
to taste-changing interactions, which is a pure lattice ar- 
tifact. Let us consider a vertex where a gluon carries a 
momentum k+Hc where lie = if / a-C and a quark/anti- 
quark has momentum (p + T\a)/{<1 + n^) respectively. 
Here we assume that |fc|, |p|, \q\ <C Tt/a. This vertex can 
be expressed as follows: 

V,(fe + n c ) ■K 7 (p + n 4 ,-< ? -n B ,fc + n c ) = 



{-ig)T J s{ P - q + k){ lvd ® z d ) AtB 
i + (-i) c - 1 

^\ '- +0(ka) h^(k + U c ) 



(58) 



Here, obviously we need to choose C„ = 0. In other 
words, the longitudinal mode is not allowed to carry high 
momentum in the gluon vertex mainly because this is 
unphysical and violates helicity conservation. In the case 
of unimproved (naive) staggered fcrmions, h^ u (k+Hc) = 
&ii,v regardless of He. Hence, it is permissible to mix with 
the taste (we call this off-diagonal mixing below) and 
the mixing coefficient is substantial 0]. m contrast, in 
the case of the improved staggered fcrmions using the fat 
links of our interest such as Fat7, Fat7 and HYP (II), 



h^(k + n c ) =O + 0(feV) 



(59) 



IV. CONCLUSION 

We have studied penguin diagrams for various im- 
proved staggered fermions at the one-loop level. The di- 
agonal mixing occurs only when the original operator has 
the spin and taste structure of (7^ £§) 1) regardless of that 
of the spectator bilinear. The main result summarized 
in Theorem 1 is that the diagonal mixing coefficient 
is identical between the unimproved staggered operators 
and the improved staggered operators constructed using 
fat links such as Fat7, Fat7+Lepage, Fat 7, HYP (I) and 
HYP (II). This is a direct consequence of the fact that 
the contribution from the improvement changes only the 
mixing with higher dimension operators and off-diagonal 
operators, which are unphysical. However, Theorem 1 
has such a limitation that it does not apply directly to 
the case of the AsqTad staggered formulation, in which 
case there is an ambiguity of choosing the fat links for the 
operators. In addition, the mixing with off-diagonal op- 
erators vanishes for Fat7, Fat7 and HYP (II). In the case 
of Fat7+Lepage and HYP (I), the off-diagonal mixing is 
significantly suppressed by the factor of [/i M(1 ] 2 . 

The results of this paper, combined with those of the 
current-current diagrams in , provide a complete set of 
matching for e'/e with all corrections of 0(g 2 ) included. 
In our numerical study of the CP violation, we adopt Fat7 
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and HYP (II) . It turns out that this choice has one addi- 
tional advantage of the absence of off-diagonal mixing in 
penguin diagrams as well as those advantages presented 
in II 0,0 
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APPENDIX A: DERIVATION OF EQ. ll35ll 

First, we define the common factor Y as 

Y = (^)ri d ^^(-i) c '- A '(^¥^) A , B ,(-i) D '- B ' 

A'B' 

^E(^®^Wi1 ( A1 ) 



AB 



Using this notation, we can simplify the identities as fol- 
lows. 



]T2sm(^).i^>(fc) = Y 

A" 

£2sin(^)-i#W(fc) = Y 

M 

^2sin(^). J ff^)(fc) = Y 

A* 



e ik-B _ e ik-A' 



pik-B' _ gik-A 



(A2) 
(A3) 
(A4) 
(A5) 



Using these identities, it is easy to derive Eq. 



APPENDIX B: FINITE INTEGRALS 

We use the following abbreviation to represent the in- 
tegration measure and its normalization factor. 



(Bl) 



Using this notation, we can express I a , lb and I c as fol- 
lows. 



(B2) 



I a = / J F 2 (p)cos 2 (p p )cos 2 (p a )sin 2 (p l/ ) 
= +11.2293(3) 



F 2 (p) cos 2 (p M )cos 2 (p«) 



16( - ln(4™V) - lE + F QO oo) ~ 40.7773(6) 
+0(m 2 a 2 ) (B3) 



F 2 (p) 2-sin 2 ( PAl )-sin 2 (p Q ) 



sin 2 ^) sin 2 (p Q ) 



= Y ( - hi(4m 2 a 2 ) - lE + F OQO o) - 9.5147(1) 



16 
"3 

+C(m 2 a 2 ) 



(B4) 



where fi =/= a v and F(p) is defined as 

1 



F(p) = 



E u sin 2 (^) + (ma)2 



(B5) 



These integrals are also given in [1J| and the results are 
consistent with each other. 
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